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A COMPUTATIONAL PROCEDURE TO ANALYZE METAL MATRIX 


LAMINATES WITH NONLINEAR LAMINATION RESIDUAL STRAINS 
By C. C. Chamis* and T. L. Sullivan* 

Lewis Research Center 
National Aeronautics and 
Space Administration 
Cleveland, Ohio 

ABSTRACT 

An approximate computational procedure is described for the analysis 
of angleplied laminates with residual nonlinear strains. The procedure con- 
sists of a combination of linear composite mechanics and incremental linear 
laminate theory. The procedure accounts for initial nonlinear strains, 
unloading, and in-situ matrix orthrotropic nonlinear behavior. The results 
obtained in applying the procedure to boron/aluminum angleplied laminates 
show that this is a convenient means to accurately predict the initial 
tangent properties of angleplied laminates in which the matrix has been 
strained nonlinearly by the lamination residual stresses. The procedure 
predicted initial tangent properties results which were in good agreement with 
measured data obtained from boron/aluminum angleplied laminates. 

KEY WORDS: Fiber composites, boron/aluminum composites, nonlinear analysis, 

approximate computational procedure, nonlinear residual strains, 
initial properties, computer program. 


*Aerospace Research Engineer 


i 




1 


INTRODUCTION 


Determining the reliability of structural components fabricated 
from fiber composites requires use of various deterministic mathematical 
models. These models, for example, may be for relating stresses to 
applied forces, stress intensities at the tips of cracks to nominal 
stresses in the component, buckling resistance to applied force, vibra- 
tion response to excitation sources, etc. Deterministic models of the 
cases just mentioned require initial tangent and strain dependent stress- 

(l)l 

strain relationships. Experimental data indicate that the presence 

of residual stress has significant influence on the initial tangent 
properties of boron/aluminum angleplied laminates. Analysis of experimental 
data shows that the lamination residual stresses (the fabrication process 
induces thermal strains) may be of sufficiently high magnitude to strain the 
aluminum matrix nonlinear ly in certain ply orientation configurations. 

When the matrix is strained nonlinearly, the stress-strain relationships 
of the la mi nate become load-path dependent. Though nonlinear response of 
composite behavior has received some attention ^ the lamination residual 
strain nonlinearity facet, including its effects on subsequent loading, has 
not been examined. It is the purpose of this investigation to describe 
a computational procedure for predicting initial tangent stress-strain 
relationships of angleplied laminates in which the matrix has been strained 
nonlinearly by the residual stress. 


^The numbers in the parenthesis refers to list of references appended to 


this paper. 
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BACKGROUND AND DESCRIPTION OF PROCEDURE 


The background leading to the procedure, its theoretical basis, and 
the computational procedure evolved are described in this section. 

Background 

The computational procedure described herein evolved from the analysis 
of test data from boron/aluminum angleplied laminates ’ . These lami- 

nates had various ply configurations and were loaded at various angles to 
their material axis of symmetry. Application of linear laminate analysis 
to these laminates predicted results for the initial tangent modulus 
(modulus of elasticity) which were considerably higher than the corres- 
ponding measured values. This observation suggested that the lamination 
residual stresses are present and of sufficient magnitude to strain the 
matrix nonlinear ly in some plies. Further analysis of these data 
indicated that the residual stress may strain the matrix at considerably 
different nonlinear strain magnitudes along the fiber direction, trans- 
verse to it, and in intralaminar shear. It was concluded, therefore, 
that a computational procedure should allow for orthotropic nonlinear 
behavior of the in-situ matrix. 

A computational procedure using the finite element method ^ is 
a logical approach for such analyses. However, this method consumes 
large computational times in nonlinear analyses. Therefore, an alternate 
approximate method was pursued. In addition, this alternate method 
could be readily implemented through the use of a presently available com- 
puter code with the addition of three equations to be described sub- 


sequently. 
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Theoretical Basis 

The computational procedure evolved is an approximate analysis* It 
consists of linear composite mechanics and a piece-wise linear laminate 
analysis to handle the nonlinear response of the in- situ matrix* Tangent 
properties at the current cumulative matrix strain level are used as 
inputs to compute the strains due to the next load increment. The equa- 
tions of force equilibrium are satisfied at the micromechanics, macro- 
mechanics, and laminate levels. 

The macro- and micro-residual strains are obtained using linear 
composite mechanics and average matrix thermal and mechanical properties. 
These average properties are used as inputs in the micromechanics level 
to generate the ply properties. The ply properties and the temperature 
difference between processing and use- temperatures are used as inputs in 
the laminate analysis level to compute the residual strains and stresses 
in the plies. The maximum residual strains in the matrix are computed 
using the ply strains and the strain magnification factors . The average 
residual strains in the matrix are computed using the following equations. 

Along the fiber direction 

Transverse to the fiber direction 


^ 7 * 22 R . “ 2/^2 (-^122 “ 122 ) I ^ 


( 2 ) 
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Intralaminar shear 

^<7rti2R " ^ 

The notation in equations (l), (2), and (3) is as follows: £. 

denotes strain; LT is the difference in temperature between processing 
and use-temperature; oL is the thermal coefficient of expansion, and ^ 
is the matrix magnification factor. The subscript VI denotes matrix, 

R residual, and j 2 ply. The numerical subscript 1 denotes direction 
parallel to the ply fibers and 2 normal to them. 

Equations (2) and (3) were obtained by assuming that the transverse 
and intralaminar shear strains in the matrix decrease from their maximum 
values as a cosine function. Determining the integrated average of this 
cosine variation of strain yields the desired result. 

The suitability of approximating the strain distribution in the matrix 
assuming a cosine variation was compared with the corresponding distri- 
bution using a two dimensional second order finite element analysis. The 
results are shown in figure 1. As can be seen in this figure, the areas 
under the two curves are almost the same. Therefore, the integrated 
average of the cosine distribution for the average matrix strain is a 
reasonable approximation. 

The average residual strains in the matrix in the various plies are 
used as inputs in the micromechanics level to compute the ply tangent 
properties with initial strains. The ply tangent properties and a small 
load increment are used as inputs in the laminate analysis level. The 
laminate analysis yields the initial tangent composite properties. 



5 


Computational Procedure 

The computational procedure is in the form of a computer program. 

In the present investigation, the computer program was generated by 
modifying an available computer program Briefly, this program 

consists of a collection functional modules to carry out composite 
micromechanics, macromechanics, and laminate analysis. The inputs to 
this computer program are constituent material properties, and composite 
geometry. The program allows the in- situ matrix to be orthotropic. 

In the modified version, the temperature dependent thermal and 
mechanical properties of the matrix are read in the program in the form 
of tables. The fiber properties are read as in reference 11. A quali- 
tative flow chart of the resulting computer program is shown in figure 2. 
Note the temperature difference is read in. Note also the provision to 
check the positive definite condition of the array of the stress-strain 
relations of the aluminum matrix. Violation of this condition leads to 
predictions of negative moduli of elasticity for both the ply and the 
composite. 

The last block in the flow chart, figure 2, is enclosed by an inter- 
rupted line. This is to note that the procedure can be readily extended 

cLu-t 

to continue the nonlinear analysis/.to mechanical and/or additional 


thermal loads. 
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APPLICATION OF COMPUTATIONAL FROCEDURE 

The computational procedure described previously and shown in the 
qualitative flow chart, figure 2, was used to predict initial tangent 
properties of boron/aluminum composites. In this section, the computed 
properties are compared with the measured data reported in references 1 
and 7. The constituent materials for these composites were 4-mil diameter 
boron fiber and 6061-0 aluminum alloy matrix. The fiber volume ratio was 
about 0.50. The laminates were made by hot pressing at 970°F. 

Input Data and Test Specimens 

The input data required for the computational procedure consisted of 
tabular data of matrix thermal properties and stress-strain data due to 
mechanical load. The thermal data of the matrix was obtained from reference 
13 and are shown in table 1. 

The temperature difference used in the residual strain computations 
in this investigation was 700 °f. The reason for this is that the matrix 
starts supporting stress of greater than 1000 psi between 800 and 700 F. 

See table 1. 

The mechanical load stress-strain data were obtained from axial tension 
and from torsion tests on thin tubes made from 6061-0 aluminum. The tube 
test section was 2 inches inside diameter, 10 inches long, and .060 
inch in wall thickness. The tube was machined to the test section dimen- 
sions from commercial tube stock allowing 1 inch at each end for gripping. 
The tubes were instrumented at their midlength with strain gage rosettes. 

The tubes were loaded in a multiaxial testing machine under tension and 
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torsion loads. Load and strain data were recorded at regular load 

intervals. These data were reduced to engineering stress-strain data 

(l4) 

using a strain-gage data-reduction program . Tensile load data 
suitable for input into the computational procedure are shown in table 2. 
The corresponding shear data are shown in table 3- 

Additional input data were boron fiber material properties, composite 
geometry, and the correlation coefficients described in reference 11. 

g 

The fiber properties used were: modulus = 57 x 10 psi, Poisson's ratio = 

.2, and thermal coefficient of expansion = 2.8 x 10 ^ in/in/°F. 

Results and Comparisons 

As was mentioned previously, the computational procedure described 
herein was used to predict initial tangent properties for test specimens 
as described in reference 1.. 

A summary of the calculated lamination residual strains in the plies 
and the corresponding average residual strains in the matrix are shown in 
table Compared to the corresponding strains in tables 2 and 3> it is 
seen that the average residual strains in the matrix in some plies are of 
sufficient magnitude to strain the matrix well into the nonlinear range. 

The matrix properties selected by the computational procedure to be 
used in subsequent calculations are summarized in table 5. The ortho- 
tropic nonlinear behavior of the in- situ matrix is clearly illustrated 
by the results in this table since ^ ^mit ^ or the m j or ity of the 

cases shown. 

Initial tangent modulus and Poisson’s ratio values predicted by the 
computational procedure with and without residual strains are compared with 
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measured data for several composites in table 6. As can be seen in this 
table, the predicted moduli from the "with residual strain” case are in 
good agreement with the measured data. The predicted moduli from the "no 
residual strain" case are considerably higher than the measured data. 

In some cases, the predicted values exceed the measured ones by as much 
as 35 percent. 

The comparisons for the Poisson's ratio values, however, are not good 
for some composites for the "with residual strain" case. The difficulty 
in comparing Poisson's ratio values arises from the fact that in certain 
cases there is a transverse strain gradient through the thickness of the 
laminate. The predicted value is an average value while the measured 
value is a local one, i.e., that on the material surface. 

To illustrate the aforementioned point, the following calculations 

were performed. For the (0+30) o and the (CH-45) laminates, the Poisson's 

- o — o 

ratios were computed for laminates consisting of all O-plies, all + 30 -plies, 
or all +4-5 plies. In these computations, the presence of residual strains 
was incorporated by using aluminum matrix properties for each ply from 
corresponding plies in table 5- The results from these calculations are 
compared below with measured values. 

Laminate Poisson's Ratio 


0 

predicted 

.25 

(°±30) s 

measured 

• 33 

(±3°> s 

predicted 

1.68 
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Laminate 


Poisson's Ratio 

0 

predicted 

.27 

(0+45) s 

measured 

.30 

<f*5) s 

predicted 

.99 


As can be seen from these comparisons, the measured Poisson's ratio 
is between the two predicted values but close to that for the all 0-plies 
as would be expected. To eliminate the transverse strain gradient, it is 
suggested that the following types of test specimens be used for measuring 
Poisson's ratio values in laminates with nonlinear residual stresses: 

(1) wide, flat test specimens to minimize or eliminate free-edge effects, 

(2) thin tube specimens, or. ( 3 ) specimens subjected to biaxial stress 
fields . 

Predicted initial tangent properties "with" and ’'without" residual 
strains are summarized in table 7 for the composites investigated. As 
can be seen in this table, the differences are substantial, especially 
for the shear modulus. 


■DISCUSSION AND IMPLICATIONS 

The computational procedure described previously was found to be a 
convenient tool to compute initial tangent properties in composites whose 
in-situ matrix had been strained nonlinear ly by residual strains. In 
all the cases investigated, one load increment was sufficient to determine 
the initial tangent properties. The computer CPU time per case was less 
than five seconds in the UNIVAC 1106. 
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In the cases examined in the present investigation, no ply unloading 
occurred in any of the composites. Had the unloading occurred, then two 
or more load increments would have been required to establish the initial 
tangent properties. The actual computing time per load increment is about 
one to two seconds. 

The procedure described herein has two important distinctions when 
compared with other available approaches. These are: (l) the approximate 

treatment of the interfiber matrix strain variation and ( 2 ) the accounting 
of orthotropic nonlinear behavior of the in-situ matrix. The main advan- 
tages of this approach are: (l) computer running time economy, (2) ease 

of input data preparation, (3) measured nonlinear stress -strain data for 
the ply under combined load is not required, and (k) the three-dimensional 
tangent properties for both ply and composite are computed routinely as 
a part of the analysis . 

In view of the differences in the predicted and measured data for 
Poisson’s ratio and since the measured values appear to be sensitive to 
the ply stacking sequence, it might be necessary to measure these values 
in combined stress fields. This will assure constant transverse strain 
through the laminate thickness . 

Additional experimental data are needed to assess the suitability of 
this computational procedure in predicting initial values for the elastic 
constants under combined loading. The computational procedure can handle 
combined loadings including bending and nonuniform temperature through the 


thickness . 
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Since the response of material with nonlinear initial strains "becomes 
load-path dependent, care should he used in selecting measured properties 
for a particular design. If initial tangent properties are selected, 
the measured and predicted data should be consistent with the anticipated 
load path. 

In view of the large residual strains in the matrix in some composite 
configurations, particular care should be taken when designing components 
from those composites which will be subjected to fatigue in order to avoid 
failure in a few cycles. 


SUMMARY OF RESULTS 

A convenient computational procedure has been described which can be 
used to accurately predict the initial tangent properties of metal matrix 
composites in which the in situ matrix has been strained nonlinear ly by 
the lamination residual strains . 

Application of this procedure to various angleplied laminates made 
from 4-mil diameter boron/606l-0 aluminum alloy composites predicted 
results which were in good agreement with measured data for the initial 
tangent modulus. This agreement was possible because the in-situ matrix 
was considered to respond nonlinear ly in an orthotropic manner. 

The results showed that predicted values for initial tangent modulus, 
not accounting for nonlinear residual strains, are higher than measured 
data by as much as 35 percent for some laminates. 
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TABLE 1. - MATRIX PROPERTIES VS TEMPERATURE 
6061-0 ALUMINUM ALLOY 


TEMPERATURE 

°F 

MODULUS, 
10 6 PSI 

POISSON'S 

RATIO 

THER COEFF EXP, 
KT 6 IN./IN./°F 

YIELD 

STRESS, 

PSI 

0 

11.3 

0.310 

12.1 

8500 

100 

11.0 

.330 

12.8 

8500 

200 

10.7 

.330 

13.5 

8500 

300 

10.4 

.345 

14.1 

8000 

400 

10.1 

.335 

14.7 

7000 

500 

9.9 

.330 

15.2 

4000 

600 

9.8 

.330 

15.8 

3000 

700 

9.3 

.330 

16.5 

2000 

800 

8.8 

.345 

17.2 

1000 

900 

8.0 

.350 

17.9 

0 

1000 

7.0 

.350 

18.6 

0 


CS-69455 


TABLE 2. ■; MATRIX STRESS/STRAIN PROPERTIES 
6061-0 ALUMINUM ALLOY 


STRAIN, 

% 

STRESS, 

PSI 

MODULUS, 

PSI 

POISSON'S RATIO 

0 

0 

11. OxlO 6 

0.366 

.024 

2 590 

10.5 

.336 

.050 

5 260 

9.8 

.334 

.078 

7 920 

8.2 

.332 

.098 

9 250 

3.6 _ 

.392 

.44 

10 500 

2. 7xl(r 

.409 

1.2 

11900 

8.8x10* 

.400 

1.5 

11 900 

9. 9x10? 

.397 

5.0 

11900 

1 OxlO 1 

.450 


CS-69453 


TABLE 3. - MATRIX SHEAR STRESS-STRAIN PROPERTIES 
6061-0 ALUMINUM ALLOY 


STRAIN, 

% 

STRESS, 

PSI 

MODULUS, 
10 6 PSI 

0 

0 

3,27 

.04 

1150 

3.44 

.07 

2450 

3. 57 

.11 

3760 

3. 56 

.15 

5090 

3.32 

.19 

6430 

2.03 

.33 

7810 

.867 

.50 

8300 

.100 

1. 00 

8300 

.010 


CS-694SZ 






TABLE 4. - THERMAL STRAINS IN PLIES AND IN MATRIX 
BORON/ALUMINUM COMPOSITES; 4-MIL FIBER; 6061-0 ALUMINUM ALLOY; 
0.50 FIBER VOLUME RATIO; AT =-700° F 


COMPOSITE 

CONFIGURATION 

PLY 

THERMAL PLY STRAINS 
IN lO* 3 IN. /IN. 

INTEGRATED AVG MATRIX 
RFSIDUAL STRAINS IN 10' 3 IN./IN. 

£ U1 

e t22 

e ll2 

£ mll 

e m22 

1 

[°s] 

0 

0 

0 

0 

0 

0 

0 

loptfL 

0 

-2.48 

-a 03 

0 

-0 

0.020 

0 

L ™ J S 

+5 

-2.53 

-7.99 

-.963 

-.043 

,119 

3.02 

[02±15l 

0 

-2.50 

-7.90 

0 

-.017 

.339 

0 

L 4 J5 

+15 

-2.86 

-7. 57 

-2.72 

-.380 

1. 30 

8. 55 

[o>80]. 

0 

-2,67 

-7.48 

0 

-.190 

1.35 

0 

L l J 5 

+30 

-3.88 

-6.28 

-4.16 

-L39 

4.76 

13.1 

M, 

0 

-3.22 

-6.56 

0 

739 

3.49 

0 

L 4 J5 

+45 

-4.89 

-4.89 

-3.34 

-2.41 

8.55 

10. 5 

Mol 

0 

-4.73 

-4.73 

0 

-2.25 

9. 12 

0 

L c 

90 

-4.73 

-4.73 

0 

-2.25 

9.12 

0 


CS-69462 


TABLE 5, - SUMMARY OF MATRIX PROPERTIES 
SELECTED BY COMPUTATIONAL PROCEDURE IN PRESENCE OF 
NONLINEAR RESIDUAL STRAINS; 6061-0 ALUMINUM ALLOY 


COMPOSITE 

PLY 

MODULI IN 10 6 PSI 

POISSON'S RATIO, 
u ml2 

E mll 

E m22 

G ml2 

MATRIX ONLY 


11.0 

11.0 

3.6 

0. 34 

[°Z ±5 ]s 

0 

11.0 

1CL6 

3.3 

.35 


+5 

10.1 

8.3 

.07 

.33 

[0 2 ±15] s 

0 

11.0 

10.2 

3.3 

.35 

+ 15 

10.1 

3.3 

.04 

.36 

[O 2 ±30] s 

0 

10.6 

4.1 

3.3 

.36 


+30 

3.2 

1.3 

.01 

.40 

i 

, 0 

0 

8.4 

.6 

3.3 

.37 


+45 

2.2 

.2 

.01 

,40 

Ms 

0 

2.4 

.2 

3.3 

.40 

90 

2.4 

.2 


.40 


CS-69457 







TABU 6. - COMPARISONS OF MEASURED AND PREDICTED RESULTS 


INITIAL TANGENT ELASTIC CONSTANTS; BORON/ALUMINUM COMPOSITES; 4-MIL DIAM FIBER; 
6061-0 ALUMINUM ALLOY; 0.50 FIBER VOLUME RATIO; AT - -700° F 


COMPOSITE 

LOADING 

INITIAL TANGENT MODULUS 10 6 PSI 

INITIAL TANGENT POISSON'S RATIO 

PLY 

ANGLE 







ORIENTATION 


PREDICTED WITH 

MEASURED* 

PREDICTED WITH 

MEASURED’ 1 



NO 

RESIDUAL 


NO 

RESIDUAL 




RESIDUAL 

STRAIN 


RESIDUAL 

STRAIN 




STRAIN 



5TRAIN 



' oc 
, O 

0 

34 

34 

34 

0.24 

0.24 

0.22 


0 

34 

34 

33 

.26 

.26 

.26 

[0 2 ±15] s 

0 

-80 

33 

21 

31 

15 

30 

13 

.28 

.20 

.34 

.20 

.24 

.20 

[0 2 ±30] $ 

0 

30 

23 

23 

.33 

.72 

.33 


30 

25 

18 

16 

.33 

.06 

.09 


-22.5 

27 

21 

22 

.34 

.31 

.33 

Ms 

0 

27 

17 

18 

.34 

.69 

.30 

-37.5 

25 

14 

16 

.30 

.20 

.24 

[0 2 ±90]„ 

0 

28 

16 

17 

.20 

.01 

.10 

L L J S 

-37.5 

20 

14 

15 

.41 ; 

! 

.11 

.11 


'MEASURED VALUES WERE TAKEN AT ABOUT 10% OF THE COMPOSITE FRACTURE STRAIN. 

CS-69464 


TABLE 7. - PREDICTED INITIAL TANGENT ELASTIC CONSTANTS 
BORON/ALUMINUM COMPOSITES; 4-MIL DIAM FIBER; 6061-0 ALUMINUM ALLOY; 
0.50 FIBER VOLUME RATIO 

[Subscript x denotes property along Toad direction, y is 90° to X; c denotes 
composite property; and sn denotes shear strain-normal strain coupling] 


COMPOSITE 

LOAD 

ANGLE 

NO RESIDUAL STRAIN 

WITH RESIDUAL STRAIN AT ■ 

-700° F 

MODULI IN 10 6 PSI 

POISSON’S 

RATIO, 

v cxy 

COUPLING 

COEFF, 

v csn 

MODULI IN 10 6 PSI 

POISSON’S 

RATIO, 

v cxy 

COUPLING 

COEFF, 

v csn 

Ecxx 

tcyy 

G cxy 

^CXX 

E cyy 

G cxy 

[0.1 

0 

33.9 

21.3 

7.0 

0.25 

0 

33.9 

2L3 

7.0 

a 25 

0 

[°2»] 

0 

33.7 

2L2 

7.0 

.26 

0 

340 

22. 5 

4.5 

.27 

0 

[0 2 ±15] s 

0 

32.5 

20.9 

7.5 

.28 

0 

31.3 

18.6 

4.5 

.34 

0 


-80 

20.7 

30.9 

7.8 

.20 

.06 

14.7 

27.6 

4.5 

.20 

.13 

[0 2 ±30] 

0 

29.5 

20.4 

8.9 

.33 

0 

23.0 

6.0 

6.1 

.72 

0 

l- C -15 

30 

25.3 

21.2 

9.3 

.33 

.22 

17.5 

7.2 

3.4 

.06 

-.40 


-22.5 

26.8 

20.8 

9.1 

.34 

.21 

20.6 

5.7 

4.8 

.31 

-.29 

1 LTx 
* 
, o 

0 

27.0 

21.5 

9.3 

.34 

0 

17.1 

4.4 

6.8 

.69 

O 


-37.5 

24.7 

23.2 

9.1 

.30 

12 

14.3 

8.1 

2.4 

.20 

-.99 

[o 2 «o] s 

0 

27.7 

27.7 

6.7 

.20 

0 

15.7 1 

15.7 

6.1 

.01 

0 


-37.5 

2 a 2 

2a 2 

11.1 

.41 

14 

13.7 

13.7 

7.4 

.11 

.05 


CS-69465 




STRAIN DISTRIBUTION BETWEEN FIBERS 

BORON/ALUMINUM COMPOSITE WITH 0.50 FIBER VOLUME RATIO 
(SECOND ORDER ELEMENT, 441 NODES) 



FIGURE 1. 


FLOW CHART OF COMPUTATIONAL PROCEDURE 



*THESE ARE THE INITIAL TANGENT COMPOSITE PROPERTIES. 


CS-69463 


Figure 2. 















